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INTRODUCTION 


A computational  program  has  been  prepared  by  the  Information 
Processing  Division  to  compute  a particular  Reradiation  Function 
w(X;  j Transfer  Function  W (w,p),  Impulse  Response  w (t,P)  and 
Reradiatod  Waveform  g (t,P)  as  defined  below  in  terms  of  the 
incident  plane  wave  pulse.  This  IBM  704  Program, designated  USL 
Program  No.  0832,  is  in  Fortran  II  language  and  is  described  in 
Appendixes  A and  B.  Similar  computational  programs  are  described 
in  USD  Tech.  Memo.  %o.  2242-111-6?  and  2242-156-67. 


Reference  (a)  contains  a description  of  the  mathematical 
rusted  and  the  theory  behind  considering  reflection  as  a 
i at i on  n henomen on . 


;rarn  compute 


for  integer  vai ue 
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Copy  available  to  DCT  d">t*3  not 
permit  fully  k ;ble  reprodu  tion 


Table  1 


Card  Formats 


Contents 


9-16 
17-24 
29-32 
33-40 
4 1-40 

to.m 

52-54 


53-60 


61-63 


al 

Rp 

X, 

Xg 

c 

V 

V 

ISKP  (set  equal  to  zero  to  compute 
Re radiat i on  F 'unction ) 

JSKP  (set  equal  to  zero  to  compute 
Transfer  Function) 

KSKF  (set  equal  to  zero  to  compute 
Impulse  Response  and  Reradiated 
Have form) . 

For  long  jobs  requiring  the  use 
of  a dump  tape  at  least  one  of 
tnc  above  options  variables  should 
not  be  set  equal  to  zero 
IC3T0P  (in  reference  to  Reradiated 
Waveform  Array  (k,t)  iiSTOP  Is 
the  number  of  times  t is  incre- 
mented ■'..'her.  ’■  iias  its  maximum 
va  ,ue. 


1--  Initial  value  of  u) 

9- -6  Maximum  value  for  u) 

7-24  Initia  valve  of  t 

■r3-32  Initial  vo.  ue  o’  > (if  not  corrputeu) 

KK  (if  set  equal  to  ?«ro,  initial 
va.. ue  of  !;  v 1 ue  computed) 
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Card  No. 

Cols. 

Content: 

2 

39-46 

Maximum  value  for  k 

3 

1-8 

y1  "V  components  of  A 

9-16 

17-24 

25-32 

>>!  \ ' 

components  of  N 

33-40 

Axt 

4.1-48 

Ax2 

49-50 

A t 

57-64 

Increment  of  u) 

65-72 

A k 

4 

1-8 

9-16 

Au> 

17-24 

Maximum  value  of 

25-32 

Initial  value  of  'b 

33-40 

AH 

41-48 

f 

Formats : 

Card  No.  1 

- Format  6F8.3, 

513 

2 

- Format  4f8.3> 

2X,  12,  2X,  F8.3 

3 - Format  9F8.3 

4 - Format  6f8.3 


Tape  Unit  No. 

3 

4 


Tape  Units  Required 

Tape  Identification 
Data  Input 

Values  for  Re radiation  Function, 
Transfer  Function  & Impulse 
Response. 
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T spe  "nits  Pegu 'red  (cont ' d) 


-i 


7 

3 

:o 

1,3 o 5 must  be  down  to  dump 
Ko  other  sense  switches  are  used 


Celeomp  Plotter  containing  values 
for  Reradiation  Function 
Rcradiated  Waveform  Array  (k,t) 
Transfer  Function  Array  (k,w) 
Impulse  Response  Array  (k,t) 

Dump  Tape 


j 


Subroutines  Required 

Subroutine  AMP  computes  the  values  for  A(t)  array  referred  to  under 
equation  (1). 


Subroutine  Besger*  computes  the  values  o:'  cylindrical  Bessel  Functions 
(See  reference  (b)  an]  Appendix  C). 


PROGRAM  OUTPUT 

Tape  Tii  contains : 

(1)  The  values  for  the  A array  plus  the  corresponding  va  ues 
for  the  Re radiation  Function  according  to  Format  (lX,  F10.5,  5X 
F10.5) 

i 

(2)  The  values  for  th  product  of  k and  v plus  the  corresponding 
values  o the  Transfer  Function  according  to  Format  ( IX,  F10.5,  5X> 
F10.5). 

k (3)  The  values  for  t/k  plus  the  corresponding  values  for  the 

Impulse  Responst according  to  Format  (1X»  Fi0.5>  5X>  F'0.5) 

% 

\ 

i Tape  "5  contain:;: 

i 

The  c foe  the  Re  radiation  Function  ( Cal  comp  Plotter  tape) 
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Tape  7,6  contains: 

The  Reradiatcd  Waveform  Array  (k,t)  according  to  Format  (F10.B) 
Tape  -//7  contains: 

The  Transfer  Function  Array  (k,w)  v.rith  Format"  (F10.5) 

Tape  #8  contains : 

The  Impulse  Response  Array  according  to  Format  (F10.5) 

Tape  $0  is  a dump  tape. 

Notes : This  program  contai ns  options  to  compute  or  not  to  compute 
any  of  the  functions  mentioned  above.  Tape  Unit  Nos.  6,7,  and  8 can 
be  used  as  input  to  USL  Program  No.  0809,  "Representation  of  Surfaces 
A Computer  Program  to  Plot  Contours  and  Draw  Perspective  Views",  by 
Edward  Beardsworth,  Jr. 


SUMMARY 

An  IBM  7 Oh  Fortran  program,  USL  Program  No.  OO32,  has  been 
written  to  compute  a particular  Reradiation  Function,  Transfer 
Function,  Impulse  Response,  and  Rcraaiated  Waveform  in  terms  of 
the  incident  plane  wave  pulse. 


J0.  CL. 

D.  A.  STREMSKY  « 
Mathematician 
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APPENDIX  A 

NOMENCLATURE  LISTING  FOR  USL  PROGRAM  NO.  O832 


S(I) 

3(1) 

yiriT  I 

RERAD(I) 

Element  of  Reradiation  Function 
Array 

TRFER(I) 

Element  of  Transfer  Function  Array 

AKW(LM,I) 

RESP(LM,l) 

Element  of  Impulse  Response  Array 

RATIO  (LM,I) 

t/k 

GSLTM  (LM,J) 

Element  of  Reradiated  Waveform  Array 

A1 

al 

A 2 

«2 

XI 

Xl 

X2 

x2 

c , 

C 

V 

V 

N 

W 

vr 

WMAX 

Maximum  value  for  v 

T 

Initial  value  for  t 

AK 

k 

P 

A 
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AKMAX 

* 

B1 

B2 

B3 

b4 

35 

B6 

B7 

i 

BS 

BIO 

312 

OMEGA 

1 

DELTA 

TT 

TAV 

PHI 

* 


Maximum  value  for  k 


components  of  N 


Ax2 


Increment  of  w 


A w 

Maximum  value  of 
Initial  value  of  f 
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APPENDIX  B 

FLOW  CHART  FOR  USL  PROGRAM  NO.  0832 
"RERADIATION  FUNCTION,  TRANSFER  FUNCTION,  IMPULSE  RESPONSE  (CASE  1A)" 
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PFBARIATIOM  FUNCTION. TRANSFER  FLECTION  * IMPULSE  RESPONSE  (CASE  1A) 

0 . A , STR£M  SK  V 

DIMENSION  Z(500) .RFRAC (500> .APPLE (500) »A«W (50.50) «B(200) *TRFER(50» 
150) *RESP( 50.50) .RAT  10 ( 50 .50 ) .5(500) *GSUM (50.50) »R(1000) .IDUMP(lB) 
DIMENSION  bUFFER(l02A) .XAXISI5C0) .VAX!S(500) 

STOJN  ALF  #0832 
STECD  ALF  « 

WRITE  OUTPUT  TAPE  4.5501 

9501  FCRMAT(IHI) 

REAC  INPUT  TAPE  3.9502.DI 

9502  format (as ) 

IF (TDjN-DI) 9503.9504.9503 

9503  PAUSF  6 

9504  WRITE  OUTPUT  TAPE  4. 9502.01 
WRITE  OUTPUT  TAPE  4.9505 

9505  FCRMAT(10X32HD.A.STREMSKY*RCCM  3124. CODE  2242) 

REAC  INPUT  TAPE  3 . 100 . A1 ,A2 .X 1 .X2 .C .V .N . I 5KP . JSKP .KSKP .NSTOP 

100  FCRMAT(6F8.3,5I3) 

Read  input  tape  3.101  .w.wmax.t.ak.ick.akmax 

101  FCRMAT(4F8,3»2X.I?.2X«F8.3> 

READ  INPUT  TAPE  3 . 103 .B1 .B2 .P3 .e4 .05 .06 .B7 »B8  .B12 

103  FORMAT (9F8, 3) 

READ  INPUT  TAPE  3.104. AMEGA.CElTA.TT.TAu.PlO.PHI 

104  FORMAT (6F8»3) 

W 1 =w 

AX  1 bAK 
Tl»T 

TAUlxTAU 
NSTCP=NSTCP>1 
/ NPlsN.l 
NP2*N*2 
C1*N*1 
PIE«3.1A15 
DEG*1B0./PIE 
C2*2.««C1 
C3*N *2 

NPRCDxI  . 

150  DC  160  I»i.NPl 
NEWI«I 

NPPCD^NPRCD^NEWI 
160  CONTINUE 
PRCDaNPROD 
C4«C2*PR0D 
IF(ISKP) 116. 102. 116 

102  1*1 

105  S(I)s(Xl/Ai)**2.(X2/A2)**2 
SX*S  ( I ) 

2(I)*S0RTF(SX) 

ZM1=S(I)-1.0 
IF (ZM1) 110.110.108 
108  ReRAD(I)*0. 

M = I-1 
GO  TO  261 

110  ApPLE<!)*l.-SU) 

I F < N ) 1 12.111.112 

111  PEAR*1.0 
GC  TC  113 

112  PFAR*APPLE ( I ) *«N 


113  «ERAC(I)*PEAH#C1/PIE 

115  X 1 sX 1 *B5 
X2*X2*B6 
I » I ♦ 1 

GC  TO  105 

261  WRjTt  OUTPUT  TAPE  4*251 

251  FCRPATUX36HZ  RERADI AT  ION  FUNCTION) 

WRITt  OUTPUT  TAPE  4*252*  (Zd)  »RERAD(I)  »I«l*Ni) 

252  FORMAT ( 1X«F10,5.5X*F1C.5) 

URITk  OUTPUT  TAPE  4*253 

253  FORMAT (/// ) 

116  LM« 1 

125  IF(JSXP) 195*126*195 

126  I«1 

IF(KK) 130*127.130 

127  PlaBl-C»B3/V 
P2«B2-C»B4/V 

Rp« ( A 1 *P 1 ) «#2* ( A2*P2 ) **2 
AK«SCRTF(RB)/C 
130  I«1 

263  VfRITfc  OUTPUT  TAPE  4*254 

254  FORMAT (1X34HKW  TRANSFER  FUNCTION) 

165  IF(AK.AKMAX) 170*170*195 

170  IF(W-WMAX) 180.180.190 
180  AKto  <LM  • I ) »AK*W 
IF(W) 185*182*185 
182  TRFER(LM*I)«1. 00000 
GC  TC  187 
185  BCPaA<W(LM*I) 

CALL  BESGEN (BOP.B) 

BANG»BCP**Cl 

TRFER (LM  * I ) «C4*B (NP2 ) /BANC 
187  1*1*1 
Wato*B8 
GO  TO  170 
190  N2»I«1 
N4*LN 

WRlTfe  OUTPUT  TAPE  4*255* ( < AKW (LM * I ) »TRFER (LM » I ) » 1*1 ,N2> *LM.N4.N4) 

255  FCRMAT(1X»F10.5«5X»F10.5) 

280  WRITt  .OUTPUT  TAPE  7 .273  * ( (TRFER (LM*I) * I • 1 *N2) *LN«N4*N4) 

273  FORMAT ( F L 0 « 9 1 
1*1 
WsWl 
LM»LM*1 
A<«AK*B12 
GC  TO  165 
195  lai 

ENC  FILE  7 
END  FILE  7 
LM  = 1 
A<=A<1 

IF(KSKP)250.197.250 
197  MprcObI 

DC  210  M«1*NP1 
MNEW»2N-1 

mprcDsmprod»mnlw 

210  CONTINUE 
PRCC2*mprco 


n 


C5*C4/  (PIE*PPC02> 

265  WRITt  OUTPUT  TAPE  4,257 

257  FCKNAT(1X34HT/k  IMPULSE  RESPONSE) 

211  ApSFT  *ARSF ( T ) 

212  IF(ABSFT-AK)230. 230*215 
215  RESP(L.M,I)*o. 

N3»I-1 
OC  TO  235 

230  RATIC(IM*I)=T/AK 
02=RATI0(LM,I)*#2 
G3»l.-G2 

FN*N 

EXP*FN*»5 
G4*G3**EXP 
RESPaV«I)»C5*G4 
T«T*B7 
1-1*1 
GC  TO  211 

235  N3«I-1 
N4sL^ 

WRITE  OUTPUT  TAPE  4 *258 . ( (RAT IC (L* * I) «RESP (LM ♦ I ) • 1*1  * NS) .LM«N4.N4> 

258  FCRMATUX*E10.5*5X*F10.5) 

282  WRITE  OUTPUT  TAPE  8 *275 . ( (RESP < LM . I ) • I ■ 1 *N3 ) »LM«N4 *N4) 

275  FCRVAT(F10.5) 

N3P1sN3*1 

lF«NSTCP-N3Pl)23A.23A*23l 

231  DC  232  I=N3Pl.NSTnP 
RESP(N4*I)»0*0 

232  CONTINUE 

WRITE  OUTPUT  TAPE  8 , 233  * < (RESP  <L* , I)  , J-N3PI  »NST©*»  *N4 ) 

233  FORMAT  (F10.5) 

234  IF ( AK-AKMAX >236*237*237 

236  A<»AK*B12 
LM»lM*l 
1.1 

TsTl 

GO  TO  230 

237  LPsl 

END  FILE  8 
EnD  FILE  8 
N5»0 
Ak*AK1 
I»1 

238  D1UB7/AX 
T*T1 

J«l 

Dl2sA0SF (Dll) 

239  GSLE1  (L^‘«J)  so. 

240  FPACTsTAU/TT 
El=APEGA*DELTA«FRACT/2.0 
E2s£1»TAU 

E3=E2»PHI 
E4=CCSF (E3/DEG) 

call  avp<tau*R> 

FCNsR ( I ) #E4 
TCIF«(T-TAU)/Ax 
GRAPEsABSF (TOIF) 

pluv=grape-i.o 


a 


IF  (PLUM  850*850*279 
279  E9»0.0 

GC  TC  PCI 
*50  E5«GRAPE/C12 
NE5.F5 
IA.NF5M 
IP«NE5*2 

F6«G>»APE-RATIC  C L • I A > 

E7sE6/DU 
E83 1 »0<»E7 

RSPN*F7*RE5P(LM»IB)*E8»RE5P(LM«IAJ 
E9sfCN*R5PN#B10 
*01  GSUM<LM.J)»GSUf'(LM.J)*E9 

tau«tau*bio 

i.I*i 

ElO«APSF (TAU> 

IF(TT-E10)24l.240.240 

241  GSUM(LM.J) «GSU^(LM,J) /AK 

absft«absf  m 

IF(AB5FT-AK)242.243*243 

242  T»T*B7 
J«J*1 
TAU»TAU1 

IFISENSE  SWITCH  519997.9999 

9997  00  9998  LK«l«15 
IDUMP (LK)  «*0 

9998  CONTINUE 

I DUMP  (16)  ■••6 

I DUMP ( 17) **0  ' ^ " 

IDUMP (18) »N5 
CALL  DUMP (IDUMP) 

9999  GO  TO  239 

243  IF (AK»AKMAX)244*244»249 

244  AK«AK«P12 
N3«J-1 
N4«Lm 

276  WRITE  OUTPUT  TAPE  6 .277. ( (GSUM (LM» I ) » I «1 »N3) »LM«N4*N4) 

277  format «fio»5) 

N3P1«N3*1 

IF (NSTCP-N3P1) 224.224*221 

221  DC  222  I«N3P1.NST0P 
GSUM (N4. I ) *0.00000 

222  CONTINUE 

WRlTt  OUTPUT  TAPE  6 .223 . ( (GSUM <LM . I ) , I «N3P 1 .NSTOP) .LM«N4 ,N4 ) 

223  FORMAT  (FlO. 5) 

224  N5«N5*J*NST0P-N3 
IF(AK»AKMAX)278, 278.249 

278  LM»LM*1 
TsTl 

TAU«TAU1 
1 = 1 

GC  TO  238 

249  ENC  FILE  6 
END  FILE  6 

250  Ml*Nl*l 

M2sNl.2 

283  IF(ISKP) 350.259.350 

259  CALL  PLOTS  (BUFFER ( 1024) .1024.5) 


cc  260  J»l.Nl 

XAX  I 5 ( J) =Z ( J> 

YAXlS(J) «RERAD(J) 

260  continue 

CALL  PLOT  (0.0.5.0.-3) 

CALL  SCALE  < TAX  I S .5.0 «Nl » I « 10. 0) 

CALL  SCALE  (XAX IS* 10.0 »Ni • 1 » 10.0) 

CALL  LINE  (XAX1S.YAXI5.N1. 1.1*11) 

CALL  AXIS  (0.0.0.0.20HRERADIATI0N  FUNCTI0N*20*5.0«90.0.YAXr3mi>  «V 
1 AX  15(^2) *10.0) 

CALL  AXIS  (0. 0, 0. 0.1HZ. -1*10. 0*0. 0. XAXIS(Nl) »XAXl5(M2> .10.0) 

CALL  °LOT  (0.0.0.0*999) 

330  READ  INPUT  TAPE  3. 9502. ED 
IF  (ED-TE0D)9503, 9509. 9503 

9309  WRITt  OUTPUT  TAPE  A.9511 
9311  FCRNATUHOEND) 

end  file  a 

9310  STOP  5 

END  (1*1 *0*1*1) 


20 


SUBROUTINE  BESGEN  <X»B) 
DIMENSION  B (200) 

BE  * l.OE-5 
Y*S INF (1.0) 

Z = CCSFU.O> 

DsSGRTF (2.0) 

DC  5 1*1.200 

5 em=o.o 

IF  (X-10.)  6*7.7 

6 RN*35,0/ (3.5-LOGF(X) ' 

GC  TO  P 

7 RN=1.05»X*25. 

8 N*RN 

IF  (N-198)  9.9.206 
206  8(200)  *0.0 
B ( 199 ) * BE 
207  AN  * N 

STORE  = 2.0*AN«B<1991/X-B<200) 
e < 20U ) * B ( 199) 

B ( 199)  * STOHE 
NsN-1 

IF  (N-198)  280.280.207 

9 B (N*2 ) *0.0 
B ( N ♦ 1 ) *BE 

280  MAX*N*2 

DO  10  1*1. N 

JsMAX-1+1 

AN=J-2 

B(J-2)*2.0»AN*B(J-l)/X-B(J) 

10  CONTINUE 

CALL  JO(X.BES) 

Cl»BES/B ( 1 ) 

DC  15  1*1. MAX 
B ( I ) *B ( I ) *Cl 

15  CONTINUE 

16  RETURN 

END (1*1 .0.1.1) 


